Oscillator strength distribution in C3H6 isomers studied with the 
time-dependent density functional method in the continuum 

Takashi Nakatsukasa a * and Kazuhiro Yabana b ' 

a Physics Department, Tohoku University, Sendai 980-8578, Japan 

b Institute of Physics, University of Tsukuba, Tsukuba 305-8571, Japan 

We present photoabsorption oscillator strengths for C3H6 molecules with emphasis on the difference between 
isomers, cyclopropane and propylene. We use an iterative numerical method based on the time-dependent local 
density approximation with continuum, which we have recently developed. The oscillator strengths for the two 
isomers differ at photon energies above their ionization thresholds. The magnitude and the shape of the oscillator 
strength distribution are in good agreement with recent experiments. The differences between the isomers arise 
from difference in symmetry of electronic states and different behaviors of continuum excitations. 



1. Introduction 

The photoabsorption and photoionization cross 
sections of molecules are of significant interest in 
many fields of both fundamental and applied sci- 
ences. The oscillator strength distribution char- 
acterizing the optical response is the most impor- 
tant quantity in understanding the interaction of 
photon with electrons in atoms, molecules, and 
matters. The oscillator strength distribution in 
the whole spectral region has been extensively 
studied with the advanced synchrotron radiation 
and the high resolution electron energy loss spec- 
troscopy |1|2, . In order to see how the oscillator 
strength changes with varying molecular struc- 
tures, it is useful to study isomer molecules. Since 
the isomers consist of the same kind and the same 
number of atoms, we expect similarity of the os- 
cillator strengths at high photon energies. This 
is because the molecular structure has little in- 
fluence on the excitation of inner core electrons. 
However, the valence photoabsorption may differ 
according to the difference of electronic states be- 
tween the isomers. In fact, Koizumi et al. have 
observed a prominent distinction for the cross sec- 
tions of simple hydrocarbon isomers, C3H 6 (cy- 
clopropane and propylene), at photon energies of 
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10 — 20 eV [3]. The photoabsorption and pho- 
toionization data of cyclopropane in this energy 
region were later improved by a measurement 
with metallic thin film windows [J. These works 
clearly show that the oscillator strengths have dif- 
ferent peaks and shoulders depending on the iso- 
mers in a continuous spectral region above the 
ionization potentials (IPs). 

Theoretical investigation for the isomer effect 
of C3H6 has been demanded for a long time, 
however, none has been reported so far. This is 
due to difficulties in treatment of the electronic 
continuum in a non-spherical multicenter poten- 
tial. There are several methods which are able 
to take into account correlations among valence 
electrons in the continuum [5I6I7I8I9] . Neverthe- 
less, some of them are not suitable for calculat- 
ing detailed structures of the spectra, and some 
are difficult to be applied to large molecules. We 
have recently developed an alternative theoreti- 
cal method for this purpose JUj. The method 
is based on the time-dependent local-density ap- 
proximation (TDLDA) in a grid representation of 
the three-dimensional Cartesian coordinate, and 
utilizes the Green's function to take account of 
the continuum boundary condition. We use it- 
erative methods to solve linear algebraic equa- 
tions for construction of the dynamical screen- 
ing field above the first IP. The theoretical back- 
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ground of our method is similar to the one of 
Rcf. in which the authors used a single-center 
expansion technique. However, the application 
was limited to small axially symmetric molecules, 
because of difficulties in the single-center expan- 
sion. Our method is based on direct calculation 
of the self-consistent screening potential in the 
three-dimensional grid representation, which does 
not rely on the expansion and requires no spatial 
symmetry. 

In the present Letter, we report the valence 
photoabsorption of the C3H6 isomers studied 
with the continuum TDLDA method in Ref . ^Oj > 
and would like to show the power of the method. 
Then, we give an interpretation of the continuous 
spectra and elucidate origins of the isomer effects. 

2. Theory and computational method 

Optical response of molecules is characterized 
by the oscillator strength, denoted as df/dw in 
the followings, which is given by 

J^ = _2ma; im ^ J d 3 rr„5n v (r,uj), (1) 

v— x,y ,2 

where the transition density 8n v is related to the 
Fourier component of a time-dependent exter- 
nal dipole perturbation in ^-direction, V v (ui) = 
r„, through a complex susceptibility (5n v {uj) = 

The TDLDA describes a spin-independent N- 
electron system in terms of the time-dependent 
Kohn-Sham equations. Correlations among elec- 
trons are taken into account through deforma- 
tions of the self-consistent Kohn-Sham potential. 
Linearizing the Kohn-Sham potential with re- 
spect to the transition density, we obtain 



5n v {r,u) = I <Fr xo(r,r ;w) 



Sn(r") 



8n u {v") 



(2) 



system without the correlations and is given by 

occ 

X o(r,r'; W ) = 2 J>(r) { (G(r, r'; u - 



(3) 



+ G(r,r'; ej + W )}^(r'). 



Here, </>j's are the ground-state Kohn-Sham or- 
bitals and G is the Green's function for an elec- 
tron in the static Kohn-Sham potential. In order 
to properly treat the electronic continuum, the 
outgoing boundary condition must be imposed on 
G. Construction of the Green's function is easily 
done for a rotationally invariant potential Vq(t), 
using the partial wave expansion jllj. Thus, we 
split the Kohn-Sham potential into two parts, 
a long-range spherical part, Vo(r), and a short- 
range deformed part, V(r) = Vks(r) — Vb(r). 
First, we construct the Green's function, Go, for 
the spherical potential Vo, then, G can be ob- 
tained from an identity 



G — Gq GqVG. 



(4) 



The xo( r , r ';w) is a complex susceptibility for a 



We solve Eqs. (0, ©, and (@J simultaneously 
in the uniform grid representation of the three- 
dimensional real space. Equations and Q 
are linear algebraic equations with respect to Sn 
and G, respectively, for which an iterative method 
provides an efficient algorithm for the numeri- 
cal procedure. We adopt the generalized con- 
jugate residual method for these non-hermitian 
problems. We would like to refer the reader to 
Ref. ^U] for detailed discussion of the methodol- 
ogy and the theoretical background. 

The exchange-correlation potential is a sum 
of the local density part given by Ref. [T3| . 
H^ pz ^[p], and the gradient correction of Ref. [Tl|. 
^( LB ) which will be abbreviated to LB po- 

tential. This gradient correction is constructed so 
as to reproduce the correct Coulomb asymptotic 
behavior of the potential (— e 2 /r) and to describe 
the Rydberg states. It was also pointed out that 
the LB potential is necessary to reproduce the ex- 
citation energies of high-lying bound states in the 
TDLDA 15 . In our previous work we have 
also found for simple molecules that the TDLDA 
with the LB potential reasonably accounts for res- 
onances embedded in the continuum. 
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Table 1 

Calculated eigenvalues of occupied valence orbitals in units of eV. Values in brackets indicate eigenvalues 
calculated using a different value of the parameter in /i( LB ) (/? = 0.05). See text for details. Line styles in 
Figs. [I] (b), [3(b), and|3|(b) are indicated in the third and sixth columns. We use following abbreviations: 
"S" for the solid, "Do" for the dotted, "Da" for the dashed, "LD" for the long-dashed, "DD" for the 
dot-dashed, and "T-" for the thick lines. 
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(2a") 2 
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(3<) 2 
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(14') 2 


-15.4 (-16.8) 


Da 


(la") 2 
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T-LD 


(2e') 4 
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DD 
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DD 
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-23.9 (-25.3) 


LD 
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(4a') 2 
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-22.3 
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Do 
Da 
LD 



3. Results and discussion 

3.1. Ground-state properties 

We fix the geometry of nuclei optimized for 
the ground state. This is based on a semiem- 
pirical method known as PM3 (T5|. We only 
treat valence electrons in the TDLDA calculation. 
Thus, we use the norm-conserving pseudopoten- 
tial |16j with a separable approximation 17 for 
the electron-ion potentials. The coordinate space 
is discretized in a square mesh of 0.3 A and we 
adopt all the grid points inside a sphere of 6 A 
radius. This results in a model space of 33,401 
grid points. 

First, we calculate the ground state of cyclo- 
propane and propylene by solving the Kohn-Sham 
equations with the exchange-correlation potential 
of Ai< pz ) +^ LB \ The LB potential ^ LB ' contains 
a parameter j3 and we adjust this value to 
make eigenvalues of the highest occupied molecu- 
lar orbitals (HOMO) coincide with the empirical 
vertical IPs (10.54 eV for cyclopropane [TS] and 
9.91 eV for propylene 19 ). The occupied va- 
lence orbitals in the ground state calculated with 
[3 = 0.015 are listed in Table Q] The HOMO 
eigenvalues are well reproduced for both isomers. 
Propylene has a geometry of the C s point group 
while cyclopropane has the group. Although 



these isomers possess equal number of valence 
electrons (eighteen valence electrons), the elec- 
tron configuration of cyclopropane is more degen- 
erate in energy because of the higher symmetry. 

3.2. Photoabsorption oscillator strength 

Now we calculate the photoresponse of the iso- 
mers. We use complex frequencies, id + iT/2 with 
r = 0.5 eV. The V plays a role of a smoothing pa- 
rameter to make the energy resolution finite. This 
also helps a convergence of the numerical itera- 
tion procedure JU|. In Figs.n]and|21respectively, 
we show the calculated photoabsorption oscillator 
strength, df /did, for cyclopropane and propylene 
in a frequency (photon energy) range of 8 — 50 eV. 
The calculations have been done with a frequency 
mesh of Aid = 0.25 eV for a region of 8 < u> < 20 
eV, and with Aw = 0.5 eV for the rest of fre- 
quencies. The oscillator strength distributions of 
the isomers are nearly identical at w > 22 eV. 
This energy roughly corresponds to the ioniza- 
tion energy of the lowest-lying a orbital. The 
df /did monotonically decreases as the frequency 
increases but has a large tail at high frequency. 
This behavior of the high-frequency tail in df /did 
is one of the characteristics of the electronic ex- 
citations in the continuum, which was also found 
in our previous studies of simple molecules |10) . 
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The molecular structure has a little influence on 
the electronic continuum in this energy region 
(lu > 22 eV). 

In contrast, in the frequency region below 22 
eV, different structures are observed among the 
isomers. The df/dw of propylene shows a single 
broad peak at w = 13 ~ 18 eV with small wig- 
gles. On the other hand, distinctive three peaks 
at u> — 11.8, 13.5, and 15.7 eV, are found in cy- 
clopropane. This difference exactly matches the 
experimental findings of the isomer effect (the 
thin solid line in Fig. ^ (a))- The energy posi- 
tions of calculated peaks are lower than the ex- 
perimental ones by about 1.5 eV. This is also 
true for propylene in Fig. [21 (a), in which the 
broad peak is shifted to lower energy by 1.5 eV 
compared to the experiment. We would like to 
mention that, if we treat the electrons as re- 
sponding independently to the external dipole 
field, we cannot reproduce the main feature of 
the oscillator strength distributions. We call this 
approximation "independent-particle approxima- 
tion (IPA)" , which corresponds to neglecting the 
induced screening potential, the second term in 
the bracket in Eq. In Figs. ^ ( a ) and|2(a), 
the IPA calculations are shown by dashed lines. 

The Thomas-Kuhn-Reiche (TRK) sum rule 
tells us the integrated oscillator strength /(oo) = 
18 in our calculations, since we only treat valence 
electrons in C3H6. In Tabled partial sum values 
of the oscillator strengths are listed and compared 
to the experiment. Again, the IPA calculation 
cannot account for the data, while the TDLDA in 
the continuum well agrees with the experiment. 
Calculated total sum values for 8 < u < 60 eV 
are 16.2 for cyclopropane and 15.8 for propylene. 
These values correspond to about 90 % of the 
TRK sum rule for valence electrons. 

3.3. Origin of the difference in photoab- 
sorption between the isomers 

We would like to discuss details of resonance 
peaks and an origin of the different behaviors be- 
tween the isomers. First, let us compare the IPA 
results for the two isomers (See dashed lines in 
Figs, n (a) and |21 (a)). In the energy region of 
10 — 20 eV, although the bulk structure is similar, 
cyclopropane shows a sharper main peak at 11.5 
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Figure 1. (a) Calculated (thick solid line) and ex- 
perimental (thin solid) photoabsorption oscillator 
strength distribution for a cyclopropane molecule 
as a function of photon energy. The dashed line 
indicates the IPA calculation without dynami- 
cal screening effects. The experimental data are 
taken from Ref. 0]. See text for details, (b) An 
energy region of 10 < uj < 25 eV is magnified 
and the total oscillator strength is decomposed 
into those associated with different occupied va- 
lence electrons. See Table ^ for correspondence 
between a line style and an occupied orbital. 
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Figure 2. The same as Fig. ^ but for propylene. 
The experimental data are taken from Ref. 
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Table 2 

Partially summed oscillator strengths for CsHg isomers. The experimental values are estimated from 
data in Refs. [3E1 



Energy range 


Cyclopropane Propylene 
TDLDA (IPA) Expt. TDLDA (IPA) Expt. 


10 cV < lj < 25 eV 
25 eV < uj < 35 eV 
35 eV < uj < 60 eV 


10.0 (11.1) 9.9 9.5 (10.9) 9.4 
2.9 (1.2) 3.3 2.9 (1.8) 3.1 
2.9 (2.3) 2.7 (1.8) 



eV and an additional peak at 13 eV. This may 
be due to higher-fold degeneracies in electronic 
eigenstates in cyclopropane. These peak struc- 
tures in cyclopropane remain after inclusion of 
the dynamical screening effects, while, for propy- 
lene, the strong peak at 12 eV seen in the IPA 
is diminished. Next, we shall examine this differ- 
ence in the dynamical screening effects. 

For this purpose, it is useful to calculate a par- 
tial oscillator strength |11I10| which corresponds 
to a contribution of each occupied orbital to the 
total oscillator strength. We display the partial 
df/du in the energy range of 10 — 25 eV in Figs.^ 
(b) and|21(b). One can see that the major contri- 
butions to df/du come from bound-to-continuum 
excitations of electrons near the Fermi level; the 
HOMO (3e') 4 and the second HOMO (le") 4 in 
case of cyclopropane, and the second (10a') 2 , the 
third (9a') 2 and the fourth HOMO (la") 2 for 
propylene. 

The sharp peaks at 11.8 eV and 13.5 eV in cy- 
clopropane originate from bound-to-bound tran- 
sitions of (3a'i) 2 and (la^') 2 electrons. These reso- 
nances are also seen in the IPA calculation. Then, 
the electron-electron correlation brings out coher- 
ent contributions of bound-to-continuum excita- 
tions of (3e') 4 and (le") 4 electrons. The width of 
the resonances becomes slightly larger than that 
of the IPA, because of the autoionization process. 
The peak at 15.7 eV is also produced by coherent 
excitations of (2e') 4 (bound-to-bound) and (3e') 4 
(bound-to-continuum) electrons. We only see a 
shoulder around 15.5 eV in the IPA calculation, 
however, the dynamical effect enhances the peak. 

In the case of propylene, we find several small 
peaks of bound-to-bound excitations in the par- 
tial df/dui in Fig. |3(b). However, the bound-to- 



continuum excitations, which mostly contribute 
to the broad resonance in 13—18 eV, behave 
rather independently, and do not produce coher- 
ent enhancement of those peaks. As a result, 
the small peaks in the bound-to-bound transi- 
tions are mostly smeared out in the total oscil- 
lator strength. 

We think that this difference in the contin- 
uum response could be attributed to the differ- 
ence in strength of bound-to-bound transitions. 
The df /dio of propylene shows a typical behavior 
of the dynamical screening effects. Namely, the 
oscillator strengths (and the peak at 12 eV) in 
the IPA calculation are significantly weakened by 
the induced screening held in Eq. (J2J . Conversely, 
those at energies above 16 eV are enhanced. This 
is because the real part of the dynamical polariz- 
ability changes its sign at uj w 16 eV, then the 
screening field changes into the "anti-screening 
field" at higher energies. In the case of cyclo- 
propane, the situation is slightly more compli- 
cated. Generally speaking, the real part of the 
dynamical polarizability changes its sign from 
positive to negative at bound resonances. Since 
the degeneracies of electronic orbitals are higher 
in cyclopropane, the bound-to-bound transitions 
have large oscillator strengths. Then, there ap- 
pears a reminiscence of bound resonance in the 
continuum region. The screening field suddenly 
drops down at energies corresponding to bound- 
to-bound transitions. This provides effective anti- 
screening effects to cause the peak structures in 
the bound-to-continuum transitions. 

Finally, we would like to comment on depen- 
dence of our results upon the parameter in the 
LB potential fi^ LB \ A choice of this parame- 
ter is rather arbitrary, since the value of (3 does 
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Figure 3. The same as Fig. Q] but for using a 
parameter (3 = 0.05 in the calculation. The cal- 
culation has been performed for an energy range 
of 8 - 43.5 eV with a mesh of Au = 0.5 eV. 



not change the Coulomb asymptotic behavior. In 
fact, if we adopt (3 = 0.05, the value proposed by 
the original paper |14j . we obtain better agree- 
ment to the photoabsorption spectra for both the 
isomers. In the case of using (3 = 0.05, the Kohn- 
Sham eigenvalues for occupied orbitals in cyclo- 
propane are indicated as values in brackets in Ta- 
ble ^ A binding energy of each orbital becomes 
deeper by 1.3 — 1.5 eV, though spacings between 
the orbitals almost stay constant. The calcu- 
lated oscillator strength distribution is shown in 
Fig.|21for cyclopropane. The disagreement on the 
peak positions are removed in the calculation. A 
bound peak at 10 eV is also well reproduced in 
the calculation. This peak consists of the exci- 
tations of (3e') 4 and (le") 4 electrons. These ex- 
citations have an almost identical energy when 
we use (3 = 0.05, while the excitation of (le") 4 
electrons is shifted to lower energy by 1 eV when 
using (3 = 0.015. Except for the bound peak at 
10 eV, the characteristic features of the oscilla- 
tor strength distribution are not changed, and we 
maintain the interpretation given above. 



4. Conclusion 

The oscillator strength distributions of C3H6 
isomer molecules are studied with the TDLDA 
in the continuum utilizing the three-dimensional 
Cartesian coordinate representation. The calcu- 
lation shows good agreement with experiments. 
The oscillator strength in the energy region above 
22 eV is almost identical among the isomers, how- 
ever, different peaks appear below that. This iso- 
mer effect is analysed by calculating the partial 
oscillator strength of each occupied orbital. In 
addition to the difference in properties of bound 
electronic orbitals, it turns out that bound-to- 
continuum excitations of electrons near the Fermi 
level behave differently between the isomers. The 
bound-to-bound transitions in cyclopropane pos- 
sess large strengths, and the bound-to-continuum 
transitions exhibit coherent peak structures be- 
cause of the anti-screening effects. On the other 
hand, in propylene, the bound-to-bound transi- 
tions are too weak to produce the anti-screening 
peaks for the continuum excitations. Although 
the molecular structure directly has minor influ- 
ence on the continuum, the difference in bound- 
to-bound transitions leads to variation in the dy- 
namical screening effects to change the continuum 
excitations. 
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